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Abstract 

We propose to perform fully dynamical simulations at small quark masses by reweighting in 
lO ' the quark mass. This approach avoids some of the technical difficulties associated with direct 

simulations at very small quark masses. We calculate the weight factors stochastically, using 
determinant breakup and low mode projection to reduce the statistical fluctuations. We find that 
the weight factors fluctuate only moderately on nHYP smeared dynamical Wilson-clover ensembles, 
Qh| and we could successfully reweight 16^, (1.85fm)^ volume configurations from m„ ~ 20MeV to 

r~| I mq ~ 5MeV quark masses, reaching the e— regime. We illustrate the strength of the method by 

calculating the low energy constant F from the e— regime pseudo-scalar correlator. 
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I. INTRODUCTION 



The steady progress of simulation techniques over the last decade (see e.g. 0, 0, Q, 0]) 
as well as new insights into the reasons for algorithmic failures ^ have profoundly altered 
the status of lattice QCD: With the latest generation of supercomputers essentially all p- 
regime points, including the point of physical quark masses, have become accessible to direct 
simulation. However, in the small quark mass regime the challenges are still considerable: 

• Large volumes are needed for the stability of the algorithms when Wilson fermions are 
used Idl. 



• Auto-correlation times increase dramatically towards the chiral limit 

• Statistical fluctuations of fermionic correlators become difficult to estimate since con- 
figurations with large contributions become rare as small Dirac modes are more and 
more suppressed. 

A possible solution to evade these problems is to avoid generating an ensemble with the 
fermionic weight of the small target quark mass but instead simulate a heavier quark and 
reweight to the desired ensemble. 

This approach solves all of the above mentioned issues: The algorithm is more efficient at 
the larger mass, and smaller volumes will be sufficient from the algorithmic point of view. At 
a larger quark mass the region of small Dirac eigenvalues is oversampled with respect to the 
target distribution and thus observables that receive large contributions there (e.g. pseudo- 
scalar correlators) will be better estimated. Previously, the Polynomial Hybrid Monte Carlo 
algorithm [sl has been used as an alternative way to achieve such an oversampling (oj, [13] . 

Since the variance of an observable O is again a field theoretical observable, its statistical 
error on an importance sampled ensemble depends only on the auto-correlation of O. On the 
other hand when reweighting is employed the error is not given by the quantum mechanical 
variance (O'^w) — (Ow)"^, with w being the normalized reweighting factor, but rather the 
statistical variance {O'^w'^) — (Ow)"^ which depends on the variance of the reweighting factor 
itself as well as its correlation with the observable of interest. If O and w are (strongly) 
anti-correlated, this controls and limits the statistical error on the reweighted ensemble. 

Nevertheless, if the overlap (in configuration space) of the generated and desired ensem- 
bles is insufficient, reweighting will break down due to the fluctuations of the reweighting 
factor. This limits the range of quark mass values that can be bridged by reweighting. 

In this paper we work with two degenerate flavors of Wilson type fermions, though 
generalization to other fermionic actions is straightforward. If we have an ensemble of 
configurations {Ui} generated at bare mass mi with the Dirac operator Di = D{U\ mi), 
we can reweight it to the ensemble that corresponds to bare quark mass m2 by assigning to 
each configuration a weight factor 

^ Dl[U,]D2[U^] 

Wi (X det — I , ( 1 ) 
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and calculating expectation values as 



{Oh 



E^WiO[Ui] 



(2) 



Since the reweighting factors and their fluctuations will increase with the volume, reweighting 
becomes ineflicient on very large lattices. In practice we flnd that the fluctuations can 
be controlled and quite large volumes can be reweighted. Nevertheless, reweighting is a 
technique that is most useful at small quark masses and moderate volumes - like in e-regime 
calculations. 

There are two sources of fluctuations for the weights. One is due to the small eigenmodes 
of the Dirac operator. These physical infrared eigenmodes contribute to the weight factor 

as 



The suppression of configurations due to the small eigenvalues is physical. The weight factors 
control exceptionally large contributions to quark propagators that arise on configurations 
with small eigenmodes. Reweighting in fact reduces the statistical fluctuations of many 
observables when compared to the partial quenched case. The ultraviolet, large eigenvalue 
modes are not physical, but due to their large number they can dominate the fluctuations. 
Some of these fluctuations can be removed by smoothing, and with nHYP smeared Dirac 
operators [sf reweighting is possible also in bigger volumes. However, even on an nHYP 
smeared gauge background, the UV fluctuations are still large. They are also closely corre- 
lated with the fluctuations of the nHYP smeared plaquette, giving us yet an other option 
to reduce the noise by absorbing it into the gauge action: Including a term proportional to 
the smeared plaquette in the Lagrangian reduces the UV fluctuations of the weight factors. 
This latter reduction is not essential, especially at smaller mass shifts, but extends the reach 
of the method at the expense of introducing a (very) small pure gauge term to the action. 

Calculating the determinant in Eq. ([1]) to any reasonable accuracy can be very costly. 
Fortunately it is not necessary to do that, a stochastic estimator is sufficient. In Sect. [TTl 
we describe the stochastic reweighting, and several of its improvements. Sect. IIIII describes 
the numerical tests and efficiency of the reweighting technique, and in Sect. IIVI we present 
physics results using reweighted configurations. 

II. STOCHASTIC REWEIGHTING 

When the Dirac operator corresponds to Wilson or clover fermions it has the form^ 




(3) 



D[U] = 1- kM[U] , 



(4) 



here M is a suitable combination of the hopping and clover terms 
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where k is the hopping parameter k = (2m + 8) ^. Reweighting a configuration from ki to 
K2 requires a weight factor 

w = det — J — - = det^^(n) , 

dIdi 

n = d:,^d^dI{dI)-\ (5) 

The determinant can be calculated as an expectation value 

but obtaining a reliable estimate for w is expensive, especially when Vt is not close to one. 

An alternative way is to calculate only a stochastic estimator of the true weight factor and 
do the average over the ^ fields together with the configuration average. A similar approach 



is used in the stochastic global Monte Carlo update [uJ, ll3|- We start by writing the 
expectation value of an operator 0[f/] at mass m2 as 



(0)2 = ^ (vUe-'^det{DlD2)0[U] 
Z2 J 

'det{D\Di)det-^{n)0[U] (6) 



^2 
1 



Z2 

Z2 



where 



Z2 ^ JVUe~''^det{DlD2) 

Zi JVUe-^^det{DlDi) ^ ^ 

and the expectation value is with respect to both the U and ^ fields at mi. If we consider 
the configuration ensemble {Ui, with the gauge configurations from the original sequence 
and generated independently with weight e~^^^, the expectation value becomes 

(0)2 = , (8) 

i.e. the weight factors are replaced by a single estimator 

= e-«'(^[^>]-i)«' . (9) 

The obvious advantage of the stochastic approach is that the averages of the noise sources 
and the gauge fields commute. Without introducing any systematic error, we therefore 
need only one estimator of the weight on each configuration. The disadvantage is that a 
single estimator might fiuctuate too much and introduce large statistical errors in Eq. (jH]). 
Fortunately there are several methods that can reduce the fiuctuations of Sj to acceptable 
levels. 
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A. Improving the stochastic estimator 



In this section we describe two methods, the determinant breakup |14j . [l5| and the low 
mode subtraction that we use in calculating the stochastic weight factors. Both methods 



were used in [l3| in a different context. It will be useful for the discussion to rewrite Eq. ([5]) 

as 

= ((1 - x){DlY^ + - x)D^^ + x) , (10) 

where x = Ki/k2. We consider reweighting to smaller quark masses so x < 1, though 
everything works for x > 1 as well. At leading order {1 — x) oc (mi — 1712), i.e. apart 
from a gauge field independent additive constant x"^, Q oc (mi — m2). The exponent of the 
stochastic estimator in Eq. now can be written as 

et(^] - l)e = (1 - xn^iDlr'D^'^ + x(l - x)^^iiDl)-' + D^')^ + ix' - , (11) 
which requires one inversion of the Dirac operator D2 on ^. 



1. Determinant Breakup 



The determinant breakup is based on Ref. 15(] and has been used extensively in dynamical 



simulations. The idea is to break up the interval (ki 1^2) to sections ki ^ ki + Ak — >■ 
.... ^ Ki + NAhi = K,2, and write the determinant as the product of N terms, each on a Ak 
interval. Thus the stochastic estimator takes the form 

N 

exp(-E^i(^n-l)^-) ' (12) 



n=l 



where fi„ is the analogue of Eq. ([5]) on the Hn ~^ i^n+i interval. Since the operators Qn are 
much closer to one then the original Q, the estimator in Eq. f|T2l) has reduced fluctuations. 
Calculating it on intervals requires A^ applications of D~^, so it is expensive. On the 
other hand it predicts a stochastic weight factor for the configuration at any intermediate 
Kn value, so the same calculation can be used to reweight to many different mass values. 



2. Separating the low eigenmodes 

Eqs. ([3]) and flTT]) clearly show that the low eigenmodes of the Dirac operator not only 
give large contributions to the reweighting factor but they can dominate the stochastic 
fluctuations as well. This latter problem can be reduced by explicitly calculating the low 
eigenmodes and removing them from the stochastic estimator. It is not necessary to work 
with the exact eigenmodes of Q, subtracting approximate eigenmodes works as well. 

Assume that P is an arbitrary but Hermitian projection operator, = P, P'^ = P, and 
P = 1 — P is the complementary projector. Any operator can be decomposed as 




pnp pnp \ 

PVtP PVtP J 
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1 \ / pnp o \ /i R\ 

T l)[ q)[o 1 ) 



with T = PnP{PQP)-\ R = {PQPy^PQP and 

Q = Pnp - pnp{pnp)-^pnp . (i3) 

Now the determinant can be written as 

det = det (PfiP) det Q , (14) 

so the sectors projected by P and P are separated. If the projection operator P is built 
form the eigenvectors of Q, the second term in Eq. (fT3|) vanishes and Q = PQP. If P is 
built only from approximate eigenmodes, both terms in Q are present but the second term 
gives only a small correction. We will refer to it in the following as correction term. 

Separating Dirac eigenmodes: Due to the 75 Hermiticity of the Wilson Dirac operator the 
eigenvalues of D come in complex conjugate pairs and the eigenvectors are 75 orthogonal. 
The eigenvectors of the massless operator -DqI'^a) = M"^^) the eigenvectors of the massive 
one as well, with eigenvalues A + m. The 75 orthogonality implies that one can normalize 
the eigenvectors such that (fA'|75| ^^a*) = ^a'a and the operator 

P\ = |^^a)(^^a*| 75 

is a projector. Since Px is not Hermitian, the discussion of the previous paragraph does not 
apply. Nevertheless if P\ is built form the eigenvectors of Dirac operator, Eq. fll4p is valid 
with vanishing correction term. If P projects to a few low energy modes of D2, the first term, 
det{PflP), can be calculated explicitly, and the second one, det{PflP), can be estimated 
stochastically using Eq. flTTl) . Since the eigenvectors of the massless Dirac operator can be 
used at all mass values, the overhead of the low mode subtraction is a one time calculation of 
the eigenmodes. If those modes are subtracted during the inversion, some of the additional 
cost can be recouped by the improved convergence as well. 

Separating Hermitian eigenmodes: An alternative approach is to construct the projection 
operator from the eigenmodes of the Hermitian operator 

Eqs. ([5]) and ffTOj) in terms of the Hermitian operator become 

= (l-x)2ff2-2 + x(l-x)(75i^2~' + ^2"S5)+X^ (15) 

Now the projector is Hermitian, but it is not constructed from the eigenmodes of Q so the 
correction term in Eq. (fT3!) is necessary. 

Working with n normalized eigenvectors of H2\wi) = rii\wi) the projector is P = 
J27=i \wi){wi\ , and 

PQP = ujij\wi){wj\ 
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Figure 1: The topological charge of the original configurations. 

is an n X n matrix with coefficients 

S ' 

Uij = (1 - + x(l - x){wi\-f5\wj){r]^^ + ?7j^^) + x^Sij . (16) 

Vi 

Both the determinant and the inverse of PQP is easily calculable. The leading term of the 
stochastic estimator, PQP requires the evaluation of if^^-P|^). Since the inversion is done 
on the subspace that is orthogonal to the low eigenmodes, this can be considerably faster 
than calculating H2^\^). The correction term does not require a new inversion if one uses 
the identity H2'^\wi) = l/rii\wi). 

The eigenmodes of H2 have to be recalculated for every k interval. While the initial 
calculation of the eigenmodes is usually expensive, changing k by a small amount does not 
effect the eigenmodes much, and starting from nearly correct eigenmodes the calculation 
converges fast. An alternative is to use the eigenmodes of H of the first interval throughout, 
but than one needs to calculate on every interval. The inversion on the low modes 

is expensive and we found recalculating the eigenmodes on each interval is a better option. 

III. NUMERICAL TESTS 

We have tested the reweighting on a set of 16^ configurations generated with a 2-fiavor 
nHYP clover action. At the original parameter values, (3 = 7.2, ki = 0.1278, the lattice 
spacing is a = 0.115(3) fm and the PCAC quark mass is mpcAc ~ 20MeV. We have 180 
thermalized configurations separated by 5 trajectories. The topological charge, as measured 
by an overlap operator based on nHYP Wilson fermions {Rq = 0.6 and no clover term [l6| ) 
fluctuates evenly, suggesting that the auto-correlation of these lattices is small even for this 
traditionally slowly changing quantity (Figured]). 

Figure [2] shows the Hermitian gap, i.e. the histogram of the absolute value of the lowest 
Hermitian eigenmode of the configurations 0]. Direct simulations are safe when the left 
edge of the gap is far from zero. In our case one could probably lower the quark mass to 
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Figure 2: The gap of the Hermitian Dirac operator of the original 16^, k 
configurations. The dashed Une indicates the PC AC quark mass. 



0.1278, p = 7.2 



about aruq = 0.009, but not much more. With reweighting, on the other hand, one can go 
to much lower quark masses. Configurations with near zero eigenvalues will be suppressed, 
just as they should be, so exceptional configurations do not cause problems. Of course at 
some point one has to worry about the chiral symmetry breaking lattice artifacts of such 
actions. 

We reweighted the ensemble up to K2 = 0.1281, or mpcAc ~ 5 MeV. At this value 
there is one configuration in the ensemble with negative, and a few with very small real 
(Dirac) eigenvalue. Reweighting to such a small mass is interesting not necessarily for its 
physical importance but rather to see the suppression of the determinant at work. Since the 
stochastic reweighting automatically gives the weights at intermediate mass values, we have 
done weighted spectrum calculations at k = 0.1279 and k, = 0.1280 as well. These values do 
not cut into the gap and give physically more reliable results. 

Our goal with these configurations is to probe the e-regime with Wilson fermions. Already 
the configurations at k = 0.1278 are controlled by the finite volume. From the PCAC 
quark mass and the pseudo-scalar correlator on these runs, and from our previous 16^ x 32 



dynamical runs 17|, |l8j we estimate rriTrL f» 2.8(1) on these 16"^ lattices. In the e— regime 
the eigenmodes of both the Dirac and Hermitian Dirac operators are pushed away from 
zero, while the finite volume has little effect on the PCAC quark mass. The ratio of the 
median of the gap, fl, and mpcAc increases as one approaches the e— regime. In infinite 
volume the ratio is the renormalization factor Z^, which we expect to be near 1 with nHYP 
fermions. On large but finite volume the ratio might not have any physical meaning, but 
there is indication that in practice it remains close to Za |6|]. Our large volume runs give 
fi/mpcAc ~ 0.82, while the ratio here is 1.1, showing the effect of the finite volume. (In [sl 
we quoted Ji/mpcAc ~ 0.91. Those runs on 12^ x 24 lattices with La = 1.5fm also have 
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strong finite volume effects even with the relatively heavy rriq = 70MeV quark mass.) 



A. Reweighting with Hermitian eigenmodes 

With Hermitian subtraction one has to balance the cost of calculating the low energy 
eigenmodes with the improved convergence of the conjugate gradient iteration. The cost 
of the latter is proportional to the number of intervals between the starting and ending k 
values, which in turn determines the statistical fluctuations of the stochastic estimator due 
to the ultraviolet modes. 

We have found that removing more than 6 eigenmodes did not substantially increase the 
convergence of the inversion, nor did it decrease the fluctuations of the estimator. Sub- 
tracting only 3 eigenmodes resulted in a signiflcantly more expensive inversion that quickly 
overtook the cost of calculating the eigenmodes. Therefore we have settled on subtracting 6 
Hermitian eigenmodes. Next we considered the optimal number of steps in the determinant 
breakup. This should be such that the stochastic fluctuations of the weight factors are small 
compared to the fluctuations between the weight factors of the different conflgurations. We 
found that 99 intervals between k = 0.1278 and k = 0.1281 was more than sufficient to 
achieve that. Going to smaller breakup might have worked but at some point the start-up 
cost of the eigenmode calculations dominate the cost. 

For the computation of the eigenvalues and eigenvectors of the Hermitian Dirac operator 



we use the Primme package of McCombs and Stathopoulos 19|, |20 



B. Reweighting with Dirac eigenmodes 

We have also tested reweighting with Dirac eigenmodes. Using the ARPACK package we 
calculated 20 eigenmodes and separated ~ 16 real or complex conjugate pairs. While these 
eigenmodes work on every k interval, the conjugate gradient inversion is still expensive 
(we project on D2 modes but invert the operator D2D2) and we found this approach more 
expensive than removing Hermitian eigenmodes. Our tests were done with single precision 
Dirac eigenmodes and it is quite possible that with better eigenmodes the Dirac eigenmode 
separation becomes competitive. 



C. Reweighting factors 

In this section we present results obtained using Hermitian eigenmode separation. 

Figure [3] shows the reweighting factors at three different k values starting form the original 
K = 0.1278. As expected, they fluctuate more and more as we reweight to smaller and 
smaller masses, and the last case, k = 0.1281, is just barely acceptable. At that point there 
are several conflgurations that have very small Hermitian eigenvalues, and the extreme 
suppression of those conflgurations is evident. 
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Figure 3: Reweighting factors at k = 0.1279, 0.1280 and 0.1281 on 180 16^ configurations. The 
reweighting factors are normalized such that their average is (s) = 1. 

There is another source for the large fluctuations, the ultraviolet noise. The fluctuation 
of the nHYP plaquette is a good representative of the UV noise, and it correlates closely 
with the UV part of the weight factor (top panel of Figure H]). We deflne the UV part as the 
weight factor without the explicitly calculated low eigenmodes, i.e. the value determined 
by the stochastic process. While this deflnition is not unique, it captures the correlation 
with the nHYP plaquette and suggests that at least some of the UV noise can be removed 
by introducing an nHYP plaquette term in the new action by reweighting from Si = Sg — 
lndet{DlDi) to 

S2 = Sg- In det{DlD2) + ^ E(3 - Re Tr f/p,nHYp) • (17) 

The correlation shown in Figure H] can be captured by an nHYP plaquette term with f3 = 
-0.00133 coefficient at k = 0.1281, or (3 = -0.00116 at k = 0.1280. The new term is a 
local, pure gauge term with very small coefficient. While such a term is not necessary, it 
does help in reducing the fluctuations. The lower panel of Figure H] shows the reweighting 
factors at k = 0.1281 both with (solid line) and without (dotted line) an nHYP plaquette 
term in the action. 

Of course the variance of the reweighting factor is only one aspect. The real test is how 
they combine with observables to give fully unquenched results. We present two examples 
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Figure 4: The logarithm of the UV part of the reweighting factor versus the smeared nHYP 
plaquette (top panel), and the reweighting factors after the observed correlation is removed by an 
nHYP plaquette term in the action (lower panel). The dashed red line in the lower panel is the 
original reweighting factor. The data is for k = 0.1281. 



here, both for reweighting to k = 0.1280. 

Figure [5] shows the scalar correlator for the partially quenched, reweighted and nHYP 
reweighted ensembles. The partially quenched data shows that the correlator becomes neg- 
ative, a well known "artifact" of partial quenching. This disease is cured by reweighting. 
On the fully dynamical ensembles the propagator is positive with both kinds of reweighting. 
In fact the two reweighted ensembles are hardly distinguishable, though the errors on the 
nHYP ensemble are about 25% lower. 

In our second example we look at the pseudo-scalar correlator at t = nt/2 = 8, estimated 
from point-to-point propagators at a single time slice on the individual configurations. Figure 
[6] shows the partial quenched and reweighted values, C^7r(^</2) and SjC^7r(n(/2), with weight 
factors Si corresponding to a new action 5*2 both with and without the nHYP plaquette term 
(see Eq. fll7p ). Reweighting removes the very large spike (corresponding to a configuration 
with a nearly zero eigenmode), and reweighting with the nHYP plaquette term reduces the 
fluctuations by an additional factor of two (observe the scale difference). The reweighted 
data has considerably smaller statistical fluctuations than the partially quenched one. 
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Figure 5: The scalar correlator at k = 0.1280 for the partially quenched, reweighted and nHYP 
reweighted ensembles. The reweighted data points are slightly offset for clarity. 
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Figure 6: The pseudo-scalar correlator at t = nt/2 without reweighting and with reweighting. The 
third panel shows the correlator with reweighting that includes an nHYP plaquette term. Observe 
the scale difference for the last panel. All data are for k = 0.1280. 
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K 


a mpcAC 


F(MeV) 


0.1278 


0.0119(5) 


79(3)(4) 


0.1279 


0.0090(3) 


79(4)(4) 


0.1280 


0.0062(3) 


81(8)(3) 


0.1281 


0.0027(5) 


78(10)(1) 



Table I: The PCAC quark mass and the low energy constant F on the original and reweighted 
ensembles. The first error of F is statistical, the second is systematic, due to the uncertainty of 
the parameter mTiV . 



IV. PHYSICAL RESULTS 



Our goal in this paper is to illustrate the effectiveness of the reweighting method. The 
physical results we present in this section are preliminary, they merely serve to illustrate the 
power of the reweighting technique. 

As we have mentioned in Sect. Illlt the original ensemble consists of 180 IQ'^ configurations 
generated with 2 degenerate flavors of nHYP clover fermions at coupling 13 = 7.2, k, = 0.1278. 
The lattice spacing is a = 0.115(3) as calculated from the Sommer parameter ro/a = 4.25(10) 



2ll . |22| using rg = 0.49 fm, and the PCAC quark mass is ampcAc = 0.0119(5), translating 
to mpcAC ~ 20MeV. From the quadratic time dependence of the axial correlator 2^, on 
our volume of V = (1.87fm)^ we estimate mHV ~ 2.1, in or close to the e— regime. 

Reweighting has no observable effect on the lattice spacing, both with standard and 
with nHYP reweighting ro/a = 4.25 at every k value, though the error increases to 0.15 at 
K, = 0.1281. One reason for the relatively large error is the short time extent of the lattices. 

Table|T]lists the PCAC quark mass as obtained on the reweighted ensembles. These values 
correspond to quark masses between 5 and 20MeV. mpcAc shows a linear dependence on 
the bare quark mass with Kcr = 0.1282, as shown in Figure [3 With reweighting we were 
able to decrease the quark mass by a factor of 4. 

In the e— regime NLO chiral perturbation theory predicts a quadratic form for the meson 



correlators. For the pseudoscalar meson the prediction to C(e^) is |2J, l25|, l26| 

Gpit) = rf=^x(P(x)P(0)) 



where 



ap = ^/i(2mSFp) (19) 



bp = |^(l-^/i(2mSVp)) 

are related to the chiral low energy constants S and F, while 

3A 

p = 1 H 1= 

2F^VV 
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Figure 7: The mpcAC quark mass as a function of the bare mass, mq = — l/Kcr)/2 with 
Kcr = 0.1282. 



is the shape factor with /5i = 0.14046 for our symmetric geometry 
in terms of Bessel functions, Ii{u) = 8Y' (u) / {uY (u)) . The function 



24| . Ji can be expressed 



2^^ 2^ 12^ 



(20) 



describes the quadratic time dependence. Our data follows this expected functional form in 
the region t G [4, 12]. Figure [H] shows the pseudoscalar correlator on the original as well as 
on the reweighted data sets with the corresponding quadratic fits. 

The fit gives the constant term with 6-8% error, predicting S^/'^ with 1% accuracy, 
while the ratio ap/bp has 8-30% errors, predicting F at the 4-15% level. The large errors are 
not surprising as the quadratic term measures only finite size effects, and our lattices are 
not small. The errors are particularly large at k = 0.1281 where the configuration overlap 
with the original ensemble is getting small. Nevertheless we chose to present results based 
on the ratio ap/bp as it is free of renormalization factors. To predict the low energy constant 
S we would need the renormalization factor of the pseudoscalar density. This calculation 



is in progress and we will report the results in a forthcoming publication 23| . In Table [T] 



we list the predictions for F as obtained from the pseudoscalar correlator. The first error is 
statistical, the second is systematic from the uncertainty of the quantity niLV . The central 
values correspond to mTy = 2.1 at k = 0.1278, and at the other n values we rescaled 
this according to the PCAC quark mass. It is satisfying that the values we obtain are 
consistent with each other, suggesting that all four data sets are governed by the e-regime 
predictions. The value is also consistent with recent p-regime overla p a ction calculations, 
though somewhat smaller than overlap e-regime calculations 27|, |28|, |29|. It is possible to 



determine F from the eigenvalue distribution of the Dirac operator at imaginary chemical 
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Figure 8: The pseudoscalar correlator on the original and nHYP reweighted data sets. 



potential, giving consistent results, though with large finite volume corrections 30|. Other 



meson correlators can be used in similar way to predict the low energy constants of the 
chiral Lagrangian. 

The advantage of using Wilson-clover fermions is that it is relatively inexpensive to create 
even large volume configurations. With the reweighting technique it is possible to probe a 
whole range of mass values and approach the e— regime without independent simulations. At 
present we are running simulations on 24^ volumes at the same lattice spacing (L = 2.8 fm) 
at approximately 8 MeV quark masses. Our tests indicate that reweighting to 2-3 MeV 
quarks does not introduce large statistical errors, therefore we will be able to tests the finite 
volume dependence of the low energy constants. 



V. CONCLUSION 



Dynamical simulations with light quarks are still computationally expensive and have 
to overcome technical difficulties due to large auto-correlation time, algorithmic instability 
and statistical fiuctuations. In this paper we presented an alternative to direct simulations, 
suggesting that reweighting in the quark mass to reach the desired light mass value might be a 
better alternative. We described stochastic reweighting and presented several improvements. 

Our numerical tests show that with reweighting one can easily approach the e— regime 
with Wilson-clover quarks and we presented preliminary data for the low energy chiral 
constant F. Reweighting can also provide an alternative to partial quenched studies in 
the p— regime, since the statistical fiuctuations of the reweighted dynamical lattices can 
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be considerably smaller than the partially quenched data. As an example wc showed how 
the scalar correlator, known to become negative in partial quenched simulations, becomes 
positive and follows the expected theoretical form on the reweighted (and therefore fully 
dynamical) ensembles. 

Reweighting might not be efficient in large volumes, or might be limited to small mass 
differences. Our experience indicates that this is not a problem on (1.87 fm)^ volumes 
between 20 MeV and 5 MeV quarks or (2.8 fm)^ volumes between 8 and 3 MeV quarks with 
nHYP smeared Wilson-clover fermions. 
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